High energy kink in the dispersion of a hole in an antiferromagnet 
— double-occupancy effects on electronic excitations 
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Evolution of the hole spectral function along the F — (tt, tt) cut is studied in the antiferromagnetic 
state of the Hubbard model. The kink in the calculated hole dispersion, the sharp spectral-weight 
transfer between the branches, and the drastically suppressed coherent spectral weight near k = 
(0,0), as observed recently in the high- resolution ARPES studies of cuprate antiferromagnets, are 
shown to be strongly enhanced by finite-C/ double-occupancy effects. Together with the anomalous 
spin-wave dispersion observed earlier in high-resolution neutron-scattering studies, the present study 
provides further evidence of a unified description of magnetic and electronic excitations in cuprate 
antiferromagnets in terms of the Hubbard model. 



I. INTRODUCTION 

The motion of an added hole/electron in an antiferro- 
magnet (AF) is associated with a trail of upturned spins 
and scrambling of local AF order. Theoretical studies of 
the resulting spin polaron and renormalization of quasi- 
particle properties of doped carriers due to coupling with 
antiferromagnetic spin fluctuations has attracted consid- 
erable attention in the context of spin dynamics, anoma- 
lous normal state properties, spin- fluctuation mediated 
pairing, optical conductivity etcX 

Providing important insight into the nature of 
hole/electron dynamics in an AF, extensive ARPES stud- 
ies of quasiparticle properties of doped carriers have been 
carried out recently in cuprates,? including antiferromag- 
netic insulators such as Sr2Cu02Cl2, Ca2Cu02Cl2,'^''*'^ 
and high-temperature superconductors such as 
Pb,Bi2_xSr2Cu06+5, Bi2Sr2CaCu208+5^i^>aioai 
ARPES directly probes the one-particle spectral func- 
tion A(k, Lo), so it is one of the direct probes of electronic 
excitations in solids. Over the past decade, a great 
deal of effort has been invested in further improving 
this technique. These advances offer unprecedented 
high momentum and energy resolution, thus making 
new analysis methods possible and allowing a detailed 
comprison between theory and experiment. 

Most of the theoretical effort to study quasiparticle 
properties of added hole/electron in the Cu02 plane has 
concentrated on the t— J model, with double-occupancy 
excluded completely. Earlier fits of the ARPES data 
have focussed primarily on the (tt, 0) — (0, tt) direction 
along the AF zone boundary, with quasiparticle disper- 
sion calculated in the self-consistent Born approximation 
(SCBA) for the i— i'— i" — J model, which include second- 
and third-neighbour hopping terms. ^^•^'^'^'''l^i^iii'i^'i^ 
The most recent fil4 yields parameter values t — 0.35eV, 
t' = -0.12eV, t" = O.OSeV, and J = 0.14eV. 

However, with only AF spin interactions included, the 
Heisenberg model does not provide a complete descrip- 
tion of magnetic excitations in cuprate antiferromag- 
nets. High-resolution inelastic neutron scattering stud- 
ies of the spin-wave spectrum in the cuprate antiferro- 
magnet La2Cu04 have revealed a noticeable spin- wave 



dispersion along the AF zone boundaryi^ This feature 
was explained as arising from the higher-order ferromag- 
netic spin couplings generated in the strong-coupling ex- 
pansion for the Hubbard model,^'^!^^ including the cyclic 
ring exchange interaction, and thus representing finite-J7, 
double-occupancy effects on magnetic excitations. 

Spectral function and quasiparticle properties have 
also been obtained recently for the t — t' — t" Hub- 
bard model, with the self energy due to multiple magnon 
emission/ absorption processes evaluated in the noncross- 
ing (rainbow) approximationi^ Using the same set of 
Hubbard model parameters as obtained from the spin- 
wave dispersion fitj^i the quasiparticle dispersion was 
found to be in good agreement with the ARPES data 
for Sr2Cu02Cl2,'^ thus providing a unified description 
of both magnetic and electronic excitations in cuprates. 
Strong finite-?/ double-occupancy effects due to the large 
unrenormalized dispersion term 4J7^ were also pointed 
out on the F point spectral function, which was shown 
to have a strongly reduced coherent spectral weight and 
large spectral-weight transfer to the high-energy incoher- 
ent structure. 

Recent high-resolution angle-resolved photoemission 
spectroscopy (ARPES) experiments on the square-lattice 
antiferromagnet Ca2Cu02Cl2f^ with data taken along 
the F— (tt, tt) cut, have revealed an additional high-energy 
fast-dispersing feature which merges with the known low- 
energy feature having strongly renormalized dispersion of 

0.35. A strong suppression of coherent spectral weight 
was observed as the F point is approached, with a rapid 
spectral-weight transfer from the low- to the high-energy 
feature. Observed clearly above 0.8 eV binding energy, 
the fast-dispersing high-energy feature was found to track 
the unrenormalized band dispersion. 

In this paper we suggest that these novel features ob- 
served in the electronic excitations provide further ev- 
idence of finite-?/ double-occupancy effects in cuprate 
antiferromagnets, already observed in earlier studies of 
magnetic excitations. By examining the evolution of the 
hole spectral function along the F — (tTjTt) cut, we will 
show that the large unrenormalized (classical) dispersion 
term 4 for the hole, arising naturally in the Hubbard 
model SCBA analysis, has a significant effect on all key 
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features of the hole spectral function near the F point 
— drastically reduced coherent spectral weight, spectral- 
weight transfer to the high-energy incoherent structure, 
which tracks the unrenormalized band dispersion. This 
classical dispersion term is absent in the SCBA studies 
of the t ~ J model, as doubly occupied states originally 
present in the Hubbard model are completely projected 
out. 

The hole self energy for the Hubbard model within the 
noncrossing (SCBA) approximation is briefly reviewed in 
section H. Results of our calculation for the hole spectral 
function are given in section III, showing the evolution 
along the momentum direction (7r/2,7r/2) — (0,0), along 
with comparisons with the < — J model results and with 
the recent ARPES data. Conclusions are presented in 
section IV. 



II. HOLE SELF ENERGY AND SPECTRAL 
FUNCTION 

The hole self energy in the AF state of the Hubbard 
model was obtained in the non-crossing (SCBA) approx- 
imation asi^ 



Ek(c.) ^ (k|[E(k,c.)]|k) = (a* /3:)[S(k,c.)] 



/3k 



(akMq /3k-q -I- /3kt^q Q^k-q)^ 
+ f7q - E^_^ - Sk-q(^ + ^^q) 



,(1) 



which represents hole renormalization due to multiple 
magnon emission and absorption processes. Here f2q is 
the RPA- level magnon energy, Wq, Wq are the magnon am- 
plitudes on the two sublattices, and the unrenormalized 
Hartree-Fock (HF) level hole amplitudes |k) = (ak /3k) 
and energy E'^ are given by 
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in terms of free particle energy ek = —2t{coskx + 
cosfcy) = — 4t7k, the antiferromagnetic exchange split- 
ting 2A = mil, and the sublattice magnetization to. 

In order to make contact with the t — J model result, 
we consider the analytically simple strong-coupling limit 
{U/t ^ 1), where the RPA-level magnon amplitudes and 
energy are given bj*2^ 



(v^T^-i) /2 




FIG. 1: Evolution of hole spectral function along the 
{tv/2, 7r/2) - (0, 0) cut for the Hubbard model. 
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and the hole amplitudes and energy are given by 
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with J = AP/U. Shifting the energy scale {oj + A ^ 
Lu) to bring its zero at the lower band edge, and with z 
denoting the lattice coordination number, the hole self- 
energy expression (1) reduces to 



Sk(t^) = t z 



(Wq7k-q + Wq7k)^ 



k-q 



Sk- 



_q(w + fiq) 

(5) 

In comparison with the corresponding t — J model 
result the above expression for the hole self en- 

ergy in the Hubbard model involves an additional large 
unrenormalized (classical) dispersion term 4J7^ which 
has an energy scale twice that of the magnon energy. 
The Hubbard and t — J model results therefore differ 
substantially, especially near the F point where the clas- 
sical band energy 4 « 4 J has a large contribution in 
the hole spectral function 
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III. RESULTS 

The self-consistent numerical evaluation of Y.k{uj) was 
carried out on a 52 x 52 grid in k space, and a frequency 
interval Aui ~ 0.05 for uj in the range —10 < uj < 10, 
as described earlier.— In the following, we set the energy 
scale t = 1. 

Fig. 1 shows the evolution of the hole spectral func- 
tion for the Hubbard model along the nodal direction 
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FIG. 2: Comparison of spectral functions at the V point for 
the Hubbard model (a) and the t — J model (b). 



k — (7r/2,7r/2) — (0,0). The low-energy feature forms 
the well-studied strongly renormalized narrow band, with 
the strongly peaked spectral function near (7r/2,7r/2) 
representing quasiparticles which rapidly lose spectral 
weight with decreasing k. Simultaneously, there is a 
strong spectral-weight transfer from the quasiparticle to 
a strongly incoherent high-energy feature which becomes 
dominant near the F point. In the following we will high- 
light the role of the large classical dispersion term in the 
Hubbard model on this high-energy feature. The en- 
ergy separation (~ 4t) between the sharp coherent peak 
for k = (7r/2,7r/2) and the broad incoherent peak for 
k = (0, 0) is in good agreement with the dynamical clus- 
ter quantum Monte Carlo calculations for the Hubbard 
modeli^ 

Fig. 2 shows a comparison between the Hubbard and 
t — J model results for the hole spectral function at the 
F point for different J values. Presence of the large 
unrenormalized dispersion term in the spectral function 
(6) drastically reduces the coherent quasiparticle spectral 
weight for the Hubbard model, as explained below. Fur- 
thermore, this suppression of the coherent spectral weight 
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FIG. 3: The k = (0, 0) spectral function is significantly modi- 
fied by the classical dispersion term, as shown in terms of the 
intersection of ReSk(tt;) with the lines ld -\- 4 and for the 
Hubbard and t — J models, respectively, with J = 0.5. 
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FIG. 4: Comparison of the hole spectral function along the 
(7r/2, 7r/2) — (0, 0) direction obtained from the high-resolution 
ARPES study^ (a) with the Hubbard model result (b) . 



is a robust feature for the Hubbard model, whereas the 
overall structure of the spectral function as well as the 
quasiparticle weight are fairly sensitive to the J value for 
the t — J model. 

Fig. 3, reproduced from Ref. [23], clarifies the role 
of the large unrenormalized dispersion term 4 on the 
hole spectral function (6). For k near (7r/2,7r/2), the in- 
tersection of the line a;-|-4J7^ with Rel]k(w) occurs near 
the self energy peak (Fig. 3 in Ref. [23]), resulting in the 
large hole-energy lowering and the narrow coherent band- 
width. However, for k ~ (0, 0), the large classical energy 
4J shifts the intersection to the high-energy region and 
away from the self-energy peak, resulting in the strongly 
reduced coherent quasiparticle weight. Furthermore, the 
vanishing Rel]k('^) at the intersection implies that the 
small-k states track the unrenormalized band, while the 
large ImE implies strongly incoherent motion and reflects 
the strong coupling of the hole with the multiple-magnon 
string excitations. Furthermore, as the intersection in- 
volves the third branch of the self energy (Fig. 3), the 
high-energy k = (0, 0) state represents a spin polaron 
with ~ three magnon excitations. 

Fig. 4 shows a comparison of the recent high-resolution 
ARPES data reported in ref. [5] with the Hubbard model 
result for J = 0.5. For k near (7r/2,7r/2), the spectral 
function curves in (b) have been broadened with a small 
J] (Lorentzian broadening) in order to bring the spectral 
function peaks to the same scale as in the ARPES data 
for the purpose of comparison. No broadening was re- 
quired for k near (0, 0) as Iml]k(w) is already substantial. 
The vanishing spectral weight for k — (0, 0) observed 
in photoemission experiments on cuprates finds a natu- 
ral explanation in terms of the classical dispersion due 
to double occupancy within the Hubbard model. High- 
lighting the presence of string-like excitations in the AF, 
detailed comparison of the ARPES data with the t — J 
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FIG. 5: Comparison of the ARPES intensity plot reported in 
ref.[5] (upper) and the Hubbard model result (lower), show- 
ing the kink in the hole dispersion near (7r/4, 7r/4). The en- 
ergy scale is similar for both plots. The calculated dispersion 
of both the low- and high-energy features are seen to be in 
good agreement with the dispersion obtained by following the 
ARPES peak positions of the MDCs (circles and triangles). 



model results has also been recently reported,— where 
a Gaussian broadening was introduced over the whole k 
range. 

The inability of the SCBA calculations to reproduce 
properly the finite lincwidth seen in the ARPES plots for 
k near (7r/2,7r/2) has been pointed outi^ Since magnon 
states are not exact eigenstates of the AF, a finite magnon 
damping must clearly be included, especially for the short 
wavelength, high-energy modes. 

Fig. 5 shows a comparison of the ARPES intensity plot 
(upper) with the Hubbard model result (lower) of Fig. 
4. With t ~ 0.35eV, the energy scale in the calculated 
plot (5t) is similar as in the ARPES plot. Key features 
of the ARPES intensity plot — sharp spectral weight 
transfer and the kinklike feature ("waterfall") between 
the weakly dispersing low-energy feature near (7r/2, 7r/2) 
and the rapidly dispersing high-energy feature near (0, 0) 
— are in strong resemblance with the calculated result. 
Also, the energy separation At = 1.4eV) between the 
sharp coherent peak for k = (7r/2,7r/2) and the broad 
incoherent peak for k = (0,0) is in good agreement 
with the ARPES binding energy scale. Furthermore, the 
calculated dispersion of both low- and high-energy fea- 
tures are in good agreement with the dispersion obtained 
by following the ARPES peak positions of the MDCs 
(circles and triangles). The energy separation of about 
~ Lit K, 0.4eV associated with the kink in the disper- 
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FIG. 6: The calculated F point spectral function for the t — 
t' — t" ~ J model with parameters obtained in Ref. 4 by 
fitting the quasiparticle dispersion along the MBZ boundary 
(^,0)-(0,^). 
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FIG. 7: The calculated F point spectral function for the Hub- 
bard model with parameters obtained in Refs. [21,22] by fit- 
ting both the spin-wave and quasiparticle dispersions along 
the MBZ boundary. 



sion near (7r/4, 7r/4) is also in good agreement with the 
separation between the two distinctive high energy scales 
{El = 0.38 eV and E2 = 0.8 eV) reported for the high- 
temperature superconductorsji 

Turning to realistic cuprate parameters. Fig. 6 shows 
the calculated F point spectral function for the t — t' — 
t" — J model with best-fit parameters {J /t = 0.4, t' /t = 
-1/3, t"/t = 1/4.4) obtained in Ref. 4 by fitting 
the quasiparticle dispersion along the MBZ boundary 
(tt, 0) — (0, tt). As seen, the t — t' — t" — J model result also 
yields fairly sharp coherent spectral function peak. For 
comparison. Fig. 7 shows the calculated F point spectral 
function for the Hubbard model with best-fit parame- 
ters {U/t = 7, t'/t = -0.25, t"/t = 0.15, t = 0.35eV) 
obtained by fitting both the spin-wave as well as quasi- 
particle dispersions along the MBZ boundaryj^i^ It is 
clear that the Hubbard model result yields a much better 
agreement with respect to the strongly reduced quasipar- 
ticle spectral weight observed in ARPES studies. 
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IV. CONCLUSIONS 

In the context of recent high-resolution ARPES ex- 
periments, we studied the evolution of the hole spectral 
function in the (0, 0) — (tt, tt) direction and examined the 
role of the large unrenormalized (classical) band disper- 
sion term 4 in the Hubbard model 

For k approaching (0,0), the large classical band en- 
ergy ('^ 4J) was shown to result in a strong suppression 
of the coherent quasiparticle spectral weight in the low- 
energy, strongly renormalized band, and a rapid trans- 
fer of spectral weight to a high-energy incoherent struc- 
ture which tracks the unrenormalized band dispersion. 
While for k ~ (7r/2, 7r/2), the hole propagates coher- 
ently in the AF background, with a large quantum energy 
lowering associated mainly with single magnon emission- 
absorption process, the high-energy incoherent k ^ (0, 0) 
hole state was found to correspond to a spin polaron with 
multiple (~ three) magnon (string) excitations. 

For k = (0,0), the overall structure of the hole spec- 
tral function as well as the insignificant coherent spec- 
tral weight were found to be quite robust with respect to 
the magnitude of J for the Hubbard model, whereas the 
t — J model results exhibited pronounced sensitivity. The 
calculated hole spectral function for the t — t' ~ t" — J 
model, with same parameters as obtained from the quasi- 
particle dispersion fit along the MBZ zone boundary in 



the (tt, 0) — (0, tt) direction, yielded fairly sharp coherent 
spectral function peak and the overall shape was also not 
in agreement with ARPES. 

The large classical dispersion term in the Hubbard 
model clearly provided a better description of both the 
strongly suppressed coherent spectral weight observed in 
ARPES studies and of the high-energy incoherent peak 
tracking the unrenormalized band dispersion, thus high- 
lighting the importance of finite-f/, double-occupancy ef- 
fects on electronic excitations in cuprate antiferromag- 
nets. In view of the observed spin- wave dispersion along 
the MBZ zone boundary in high-resolution neutron- 
scattering studies, and their understanding in terms of 
finite-t/, double-occupancy effects, the present study pro- 
vides further evidence of a unified description of both 
magnetic and electronic excitations in cuprate antiferro- 
magnets within an effective Hubbard model description. 

As the dominant contribution to hole/electron self en- 
ergy due to multiple magnon emission/absorption pro- 
cesses involves the short-wavelength modes which only 
require short-range AF order, extension of this study 
to finite doping is of considerable interest. The high- 
energy feature characterized by sharp spectral-weight 
transfer and kink in the dispersion has been observed 
to be a universal anomaly in recent ARPES studies on 
doped cuprates, independent of doping as well as chemi- 
cal composition<ii2ii2iiiSiii 
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